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Abstract 

Neutron spectra from the reaction K~ + d — > n + ^+n were calculated in the 
energy range E c £l_ =0 — 50 MeV using coupled channel Faddeev equations 
for the description of the KNN — irUN three-body system. The aim was 
to trace the signature of the A(1405) resonance in the spectra. We found, 
that while in the direct spectra kinematic effects mask completely the peak 
corresponding to the resonance, the deviation spectrum method |l| is able 
to eliminate kinematics and produce clear evidence of the resonance. Four 
different phenomenological KN — ttTj interactions were used in order to study 
the effect of A (1405) pole position on the neutron spectra. 

Keywords: kaon-baryon interactions, Faddeev equations, kaon induced 
q \ reactions and scattering, nuclear few-body systems 
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1. Introduction 



The A (1405) resonance plays a central role in low-energy kaon-nuclear 
physics. Being a manifestation of the assumed attraction between negative 
antikaons and nucleons its observability and properties are crucial for the 
possible existence of antikaonic nuclear clusters. Its origin and structure, 
mainly its one- or two-pole nature, the position and widths of these poles, are 
subject to vivid discussions among the representatives of different opinions, 
in particular, chiral perturbation theory versus phenomenology. For a review 
of the "state of art" see e.g. [ll, |2| and references therein. From experimental 
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point of view, the clarification of this problem is hindered by the fact, that 
this state can not be reached in two-body reactions involving stable particles, 
it can be observed only in processes having 3 or more particles in final states. 
The simplest of these is the reaction 

/• K~ + d or K~ + p + n 
K- + d ->■ (tt + E) /=0 ,i + n (1) 

\(7T + E) /=1 +p 

which has the advantage, that its dynamics can be treated exactly in the 
framework of the coupled particle-channels Faddeev approach. This formal- 
ism has been applied earlier for this 3-body system mainly for calculation 
of 3-body quasi-bound states [3|, |4J, or K~d scattering length p, [g, |7J. The 
break-up reaction ([T]) was studied in the early papers [8|, |9[ (of course, to- 
gether with values for the K~d scattering length, too) but with the main 
emphasis on the possible signal of a resonance in the AN — EiV system, close 
to the UN threshold. Probably at that time the A(1405) topic was not as 
hot as nowadays. 

The only available experimental data on neutron spectra from the reaction 



(P) are those of Tan [10J from a bubble chamber experiment. What he gets 
is a large peak near the origin with no distinguishable structure in it. It 
corresponds to our "direct" spectrum P(E n ) (see later). 



The recent papers [111 If2j and [2| are devoted to the possible observation 
of the A (1405) in the reaction (TjQ). In these papers the dynamics is treated 
in single- plus double-scattering approximation, which might be justified at 
higher incident kaon energies. Our calculation is done for low-energy kaons, 
in the range E^L =0 — 50 MeV, having in mind stopped or slowed down 
kaons. 

The paper is organized as follows: in Sect. II. we give a brief description 
of the applied formalism, Sect. III. contains the details and the input of the 
calculation, our results and their discussion are presented in Sect. IV., while 
our conclusions are in Sect. V. 

2. Formulation of the problem 

The calculation is based on the coupled-channels AGS-Faddeev treatment 
of the KNN — ttJ^N three-body system. The details of this approach have 
been already described in detail in several papers |3_|, |4j, |5j, here we shall 



recall them only briefly, mainly to introduce the notations. The operator 
AGS equations for the transition operators Un read 

U tl = (1 - SijGo 1 + YtTjGoUj!, (2) 

where i,j = 1, 2, 3 are the usual pair-spectator indices, 

Tj = V 3 + VjGoT, 

are the two-particle T-operators and Gq(z) = (z — .Ho) -1 is the free Green- 
operator. 

The configuration space Ixjy^i/j) in which these operators act, apart from 
the usual Jacobi momentum variables Xjy^ contain a discrete index V{ = 
(a, (Ti), which is a combination of the particle composition index a 

a = {1, 2, 3} = {KNiNz, 7rEi7V 2 , 7rA^iE 2 } 

and an isospin label <j{ = {hi) corresponding to pair isospin Jj of the particle 
pair i and total isospin I: 



<j; ~ 



[tjtk] 1 1 






This choice of the isospin labels corresponds to the "isospin representation", 
which is useful when isospin conserving pair interactions are used. Another 
possibility is the equivalent "charge state" or "particle" representation, char- 
acterized by the 3rd component of the particle isospins: 

ffo~Oi*,*2*,*3*} or for a = 1 <r ~ {K n 1 n 2 ,K~p 1 n2,K~n 1 ,p 2 } 

These representations can be transformed into each other with the help of 
orthogonal matrices 

Bf i = {aa i \aa i ) i,j = 0,1,2,3 (3) 

composed from 6j symbols for i,j = 1,2,3 and from Clebsch-Gordan coef- 
ficients for i = 0,j = 1,2,3. Using isospin conserving separable interactions 
of the form 

^ = ^<W:>a^i (4) 



the Ti(z) operators can be written as 

i' 



uz)=J2\&)^( z m\ (5) 



Ui V- 



with r ViV r(z) being the usual (c- number) matrix, defined as: 

(r^W)" 1 = (A^)" 1 - <&\G {z)\gfy (6) 

The matrix indices z/j in Eq.flB]) consist of the particle space label a and the 
isospin label <ii = {hi). Due to isospin conservation of our interactions, 
the coupling constant matrix A is diagonal in {hi), while for particle pairs 
% capable to change their identity (KN -H- 7r£) it has non-diagonal elements 
in the particle labels (aa'). As for the matrix elements of Go, it does not 
change particle identities, thus it is diagonal in a, and if we take averaged 
masses for particles within an isospin multiplet, it is also diagonal in pair- 
and total isospin indices (/?,/). However, if physical (unequal) masses are 
used, G will be diagonal only in "particle" representation, while in the <7j 
"isospin" representation it will acquire non-diagonal elements both in /j and 
/, proportional to the mass differences. 

The equations ([2]) for the transition operators take the form 

U a = (1 - <y Go" 1 + J2 E W^ W. I with «. I = (giPoUji 

' • ' • J 3 3 3 

(7) 
Introducing the functions Xjj.(y,-) = (X^.|$ ), where y^ is the momentum of 
the spectator particle, corresponding to the pair j and |$ ) — \fdPx) is the 
initial state with the deuteron wave function |<^) and Pk - the momentum 
of the incident kaon, we get the set of integral equations: 

xixy 3 ) = (i - SaW^o) + E E / z ^(y* y^4 z - y^K^yj 

(8) 
with the kernel 

Zw i (yi,yj) = <s&\G (z)\g%). (9) 

The size of the system fl9]) can be reduced by introducing symmetric (antisym- 
metric) combinations of X-functions, with respect to interchange of baryon 



numbering. The baryon spins do not enter explicitly in this formalism, there- 
fore the total baryon spin S remains unchanged in the process (is a conserved 
quantum number). For a given S value the total antisymmetry required by 
the Pauli principle has to be ensured by the space-isospin part. Thus for 
S = (K~pp system) we have to work with the symmetric combinations 
of X-s, while for S — 1 (our K~d system) the antisymmetric combinations 
are needed. As a result, the labeling of the unknown functions z/j = (a<jj) is 
changed to \x a = (a, a a ), where a denotes a pair of interacting particles, irre- 
spectively to which original particle composition channel they belonged and 
a a denotes the corresponding isospin values. Thus we are left with X aaa (y a ) 
and a can take the values KN, NN, £iV and 7rS (ttN is missing, since we 
neglected the ttN interaction, see next section). 

The break-up transition operator Uq\ can be expressed in terms of the 
Un-s as: 

U 01 = ~(U n + Un + U 31 ) 
and the break-up amplitude reads 

Abu = (&f\U i\$o). (10) 

For the reaction under consideration the properly antisymmetrized final state 

is 

|$/) = |x^s,yjv;0-7rs) = — ^(|x 7r £i,yjV 2 ;C'"7rS 1 ) ~ |x,r£ 2 , Y Nl 5 °V£ 2 )) 

The break-up amplitude can be expressed in terms of the X-functions as 

ABui^T,, Yn'i ^tts) = 

- W E (x*e) [t^,kn(z ~ y 2 N )X RN (y N ) + t^^(z - y 2 N )X^(y N )] (11) 

- B^g XN (uy N + fx lS )r EiViEJV (2 - |x^ E - wy Ar | 2 )X EA r(x 7rS - wy N ), 

where B% t is an isospin recoupling matrix (see Eq.flS])), u,v and w are mass 
coefficients of the transformation between Jacobi momentum sets. In Eq. (II ip 
we omitted the isospin labels, the quantities are vectors (matrices) in isospin 
space. The on-shell amplitude for a given neutron energy E n depends on 
E n ,t and the isospin labels ovi3 : 

A(E n , t, (T^s) = Abu{^ttT„ Yat; C^tte) 



with 

|yj = v^n/^iv^s; l^vrsl = y^E^N ~ E n )fi nS ; t = cos(x 7rE , y n ) 

The physically observable final state corresponds to a certain particle 
composition, not to a definite isospin state, therefore the amplitude has to 
be transformed into the <To representation, using the suitable B matrix of 
Eq.fl3j: 

A(E n ,t,a ) = J2 ( B 03) mha ^ME n ,t,a^), 

where <Jq can be {7r + X~n, 7r°E°n, 7r~E + n}. The neutron spectrum is propor- 
tional to the differential cross section 

P(E n ,t,cr ) ~ — — — - = (27r) 4 /i7rsA i JV,7rS/ix,JViv^ — \A(E n ,t,a )\ 2 

dS2 X]rE a% N dE n P K 

The inclusive neutron spectrum (when no other particles are detected) is 
given by 



P{E n )=J2f dtP(E n ,t,a ) 



(12) 



The above considerations refer to the neutrons emerging from the reaction 
K~+d — > 7r+S+n; when the energy of the incident kaon exceeds the deuteron 
binding energy, neutrons are also emitted from the reaction K~ -\-d—¥ K~ + 
p + n. Their spectrum can be deduced in a similar way to Eq.( lTT|) from the 
Xk N and Xjstn functions. The allowed energy range for neutrons from the 
first reaction is (0,E K + E d + A), while for the second it is (0,E K + E d ), 
where A is the difference of the KN and 7rS threshold energies. 

3. Details of the calculation and the input 

The main purpose of the present work is to study the possible signa- 
ture (s) of the A (1405) resonance in the neutron spectra from the reaction 
([I|). In our calculation we used the two-body interactions of [5] and p], 
which are adjusted for our three-body model. They are s-wave, separable 
isospin dependent and isospin conserving interactions with Yamagouchi type 
form-factors. 

In particular, for the two-channel KN — 7rS interaction we used two 
variants, both having a one and a two pole version for the A(1405). They 



both reproduce all available experimental data on the low-energy KN system, 
the first one is fitted to the KEK data on the kaonic hydrogen Is level shift, 
while the second one reproduces the most recent SIDDHARTA data. Their 
pole positions relevant for the A (1405) are shown in Table [1] (the negative 
real parts correspond to distances from the KN threshold). 

Table 1: Pole positions of the KN — 7rS potentials (in MeV) 



1-pole 
2-pole 



KEK 



-23.6-35.6 i (1411.0-35.6 i) 
-22.2-36.3 i (1412.4-36.3 i) 



SIDDHARTA 



-6.4-46.8 i (1428.1-46.8 i) 
-14.8-57.2 i (1419.7-57.2 i) 



The numbers in Table [T] differ slightly from those given in the original 
papers [5|, |6j . The reason is, that the above ones were calculated with aver- 
aged masses and without Coulomb interaction - as they appear in most of 
the 3-body calculations, - while the fitting to the experimental data was per- 
formed with physical masses and Coulomb interaction. Since the main aim 
of the present work is to study the appearance of subthreshold resonances of 
different type in a 3-body reaction, we kept both interactions, not only the 
more advanced one. 

The triplet NN interaction is a two-term one to account for the short 
range repulsion, with form-factors fitted to reproduce the deuteron and s- 
wave phase shifts. The S = 1 EiV interaction in the 1 = 1/2 isospin state is 
complex, since it was deduced from a two-channel UN—AN interaction, while 
for / = 3/2 it is real. The ttN interaction was neglected in our calculation 
due to its weak s-wave part. 

The total angular momentum was restricted to L = since we believe, 
that for our s-wave interactions the essential dynamics can be traced in spite 
of this limitation. Keeping the interactions s-wave, the extension to higher 
partial waves is straightforward, unlike the case of inclusion of p,d, ..-wave 
interactions. 

We considered incident kaon energies in the interval E^l =0 — 50 MeV 
[pLA B ~ o — 250 MeV/c). The system of integral equations (jHJ) in the case 
of physical masses (I — 1/2 and I = 3/2 mixed) consists of 12 equations, 
while for averaged masses - of 8. 

As a numerical method we used expansion of the unknown functions on 
a cubic spline basis, for the distribution of nodes and collocation points the 



prescription of [13J was used with a slight modification to allow nonsymmet- 
ric intervals and distributions on the two sides of the break-up singularity. 
Complete convergence of the results was achieved for ~ 20 nodes in the non- 
break-up channels, while for the break-up channels ~ 30 — 35 nodes were 
necessary. Apart from the lower dimensionality of the matrices to be in- 
verted, the use of spline expansion is especially advantageous when break-up 
amplitudes are calculated, since no interpolation of the solutions is needed. 



4. Results and discussion 

We start the presentation of our results by a "by-product": the effect of 
the physical versus averaged masses of the K° and K~ mesons on the K~d 
scattering length. The inclusion of the possibility of isospin mixing due to 
this mass difference allowed us to extend the results of p, |6J in this respect. 
Our results for the two potential versions are shown in Table [2j The results 
for averaged masses coincide with those of p, |6j , while for physical masses 
they differ by a few per cent, mainly in the real part. At present level of 
accuracy of available information - both theoretical and experimental - on 
the KN interaction and ^-nuclear clusters this difference does not seem to 
be essential. However, once it might become useful to have some numerically 
reliable information on the order of magnitude of this effect. 



Table 2: K d scattering lengths of the KN 
masses of K~ and K° (in fm). 



7rE potentials for physical and averaged 





KEK 


SIDDHARTA 




averaged 


physical 


averaged 


physical 


1-pole 
2-pole 


-1.49+0.97 i 

-1.57+1.10 i 


-1.52+0.98 i 
-1.60+1.12 i 


-1.47+1.22 i 
-1.50+1.23 i 


-1.50+1.23 i 

-1.54+1.24 i 



On the other hand, when calculating neutron spectra our interest was 
focused on qualitative signals of the A(1405) in the line shapes, therefore we 
used averaged masses, what simplified the numerical work to some extent. 

We have calculated the inclusive neutron spectra P(E n ) ( 112]) for different 
incident kaon energies, both below and above the deuteron break-up thresh- 
old. For kaon incident energies E c ^l = 1, 20, 50MeV^ the results are displayed 



in Figs. [TJ |2] and El respectively, (upper left graphsjj. 



M rt _, MeV M nI , MeV 

1433 1413 1393 1373 1353 1333 1433 1413 1393 1373 1353 1333 





E , MeV 



E ,MeV 



E' m =1 MeV 

r\ 



Figure 1: Neutron spectra for kaon incident energy E c ^l — 1 MeV . Upper left graph - 
direct neutron spectra P(E n ), eq.fTTJ]); lower left graph - deviation spectrum PDEv(E n ), 
eq. trPT)) : upper right graph - single scattering deviation spectrum P^y(E n ), eq. (TT7)) ; 
lower right graph - "original" A(1405) as 7r°E° elastic cross section. 

The overall shape of the spectra is a strong peak near the origin with no 
signal of the A(1405) resonance. The direct P(E n ) spectra are practically 



2 Since our main concern is the possible trace of the A(1405) in the line shapes of the 
calculated spectra, the arbitrary units on the Y-axes of our graphs are chosen to optimize 
visibility. 



indistinguishable for the four considered KN — 7rS potentials. For kaon 
energies above the deuteron binding energy there are two modifications: the 
neutron spectrum from the K~+d — > 7r+£+n channel show a cusp at neutron 
energies E n = E t h when the KN system is at its threshold, and additional 
neutrons show up from the K~ + d — > K~ + p + n reaction in a form of a 
structureless bump between E n = and E n = E th (on the graphs it is scaled 



M iV MeV M^, MeV 

1452 1432 1412 1392 1372 1352 1332 1452 1432 1412 1392 1372 1352 1332 



-KEK1 pole 
KEK2 pole 
SIDDHARTA 1 pole 
SIDDHARTA 2 pole 




100 120 




60 80 



E ,MeV 



E. m = 20 MeV 



Figure 2: Neutron spectra for kaon incident energies E' 
the same as in FigfT] 



cm 



20 MeV The four graphs are 



down to allow to draw it on the same plot with the other neutrons). The 
reason, why the A(1405) is not seen in these spectra is essentially kinematical: 
the neutron energy in the resonance region should exceed the incident energy 



10 



of the kaon by 10 — 20 — 30MeV, depending on the resonance position, while 
in the deuteron the neutron energy (momentum) distribution is dominated 
by the low energy part. 



M H , MeV 

1483 1463 1443 1423 1403 1383 1363 1343 



M^, MeV 

1483 1463 1443 1423 1403 1383 1363 1343 



KEK 1 pole 
KEK 2 pole 
SIDDHARTA 1 pole 
SIDDHARTA 2 pole 





100 120 140 



E , MeV 



E" n = 50 MeV 



Figure 3: Neutron spectra for kaon incident energies E C £L = 50 MeV The four graphs are 
the same as in FigfT] 

In order to eliminate this kinematical "inconvenience" Esmaili, Akaishi 
and Yamazaki (EAY) [1] propose to consider instead of P(E n ) the deviation 
spectrum: 

DEV ~ p Ie) ( } 

± nonres\ J -'n) 

In Eq. (IT51) P n onres(E n ) is a non-resonant background spectrum, containing 
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the kinematics of the reaction. Let's see, how this idea can be realized in our 
case. 

Considering the zero-order iteration of our (symmetrized) system of in- 
tegral equations fl9]), the only non- vanishing X will be the inhomogeneous 
term: 

Xkn(Vn) = {9kn\®o} = (gRNlVdPK)- (14) 

This ansatz is usually called single scattering approximation. Substituting 
eq. fTHj) into eq.f lTTj) we get the corresponding break-up amplitude: 

A S Bu g {x lvTll y N ]a lvT ) = 5 , (x 7rS )r 7rSji?A r(z-^)(5f i?Ar |v3dPft-) 

= (a^yjvKs^ly^Ptf), (15) 

which is the matrix element of the two-body T operator between the initial 
and final state. It contains two-body dynamics through the T operator and 
the kinematical input: the transformation between the Jacob i-coordinates 
and the deuteron wave function. This is basically the formula, which EAY 
used to calculate the transition amplitude from the K~d atomic state to the 
irlln continuum. As for the non-resonant amplitude they suggest to replace 
T^rn hi (USD by V^,kn- 

^fcTO^E, Vn; O-tte) = (^ttS, MK^Kivl^ A), (16) 

that is, to use the Born-approximation, which contains all the kinematics. 
Thus we have three amplitudes with the properties 

- Abu — ► three-body dynamics + three-body kinematics, 

- A s ^j — > two-body dynamics + three-body kinematics, 

- A^ n — > three-body kinematics, 
and we expect, that the DEV spectra 

P DE v(E n ) = P(E n )/P Born (E n ) ; P$™(E n ) = P sm °(E n )/P Born (E n ) (17) 

will display (reveal) three- and two-body dynamics, respectively. 

It is assumed, that the E n dependence of P nonres (E n ) = P Born (E n ) is 
basically determined by the features of the initial and final states, while the 
details of the Vkn-rt, potential (within reasonable limits) influence it only 
weakly. This expectation is important, if the deviation spectrum method is 

12 



to be applied for extracting some information on A (1405) from an experimen- 
tally measured neutron spectrum. (It is assumed, that the matrix element 
ffT6|) can be calculated in an experimental group, too.) To check this antici- 
pated mo del- independence of the method, we calculated P Born (E n ) not with 
our realistic KN — 7rS interactions, but with the simplest possible separable 
potential: 



\X,rj:\V Born \x KN ) 



X 



■^+wiy is ' RN ^ N +(P i K N ) 2 



and took 



x i=o 

A ttZ,KN 



A 



1=1 _ 

■kT,,KN 



1 , Plf = p. 



1=1 



P 



1=0 

KN 



P 



1=1 

KN 



P. 



Born 



We calculated the PDEv(E n ) deviation spectra for different (3 Born values and 
incident kaon energies. Typical results are shown in Fig. HJ The shape of 




Figure 4: Effect of ^Bom on the deviation spectra 

the Pdev does not depend on the normalization of P nonr es (see the arbitrary 
choice of the \ w % t KN- s ) an d we have normalized the Pdev curves to have 
their maxima at 1. It can be seen, that for reasonable 7r£ — KN interac- 
tion range values the spectra show practically no dependence on /3 Born , thus 
confirming the presumed (approximate) model-independence of the Pdev 
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method. According to the foresaid, the rest of our calculations were done 
with (3 Born = 3/m" 1 . 

Figs. l.,2. and 3. demonstrate our main results. Apart from the direct 
spectra P(E n ) - which show no trace of the A(1405) - we show also the 
full- and single-scattering deviation spectra, too. For comparison, we plotted 
also the (hypothetical) 7r°S° elastic cross sections as calculated from the 
KN — 7rS interaction for the corresponding 7rE energy (top scale). Since this 
cross section "feels" only the 1 = part of the KN — 7rE interaction, its 
peak is usually identified with the A(1405) resonance. Thus the similarity 
of the PoEv(E n ) and cr(7r°S°) line shapes can tell us about the reliability of 
extracting information about A(1405) from the reaction under consideration. 
Obviously, the P^v s P ec t r a show more similarity with the original A(1405) 
shape than the full PoEv{E n ) spectra, since their dynamical content is more 
or less the same. For three of the four considered potentials (KEK 1, KEK 2, 
SIDDHARTA 2) a clear resonant structure can be seen in the full deviation 
spectra, however, the shapes and positions can significantly differ from their 
"originals" . The question of quantitative derivation of the A(1405) resonance 
parameters from the (measured or calculated) deviation spectra is a separate 
and probably non-trivial problem and will not be discussed here. 

As for the fourth potential, SIDDHARTA 1, its deviation spectra do not 
show any signature of its original A(1405), although there are clean resonance 
peaks in the corresponding P^ 9 V and cr(7r°S°) curves. The reason might be 
the extreme closeness of the pole to the KN threshold combined with its 
large width. 

Finally, we asked the question, under which conditions could the reso- 
nance be observed in the direct, P(E n ) spectra. For this purpose we modified 
two of the KN — ttY, interaction parameters of one of our potentials (KEK 1), 
^Rn kn an d ^■kvkni * n sucn a wa ^' ^ a t ^ e P os ition of the A(1405) remained 
at its original place, while its width could be made smaller. The results are 
shown in Figj5j 

It is seen, that in order to show up in the direct P(E n ) spectra the width 
of the original resonance should not exceed 10 — 15 MeV, while the deviation 
spectra reproduce the original resonance shape in an acceptable way. Thus 
the real A (1405) peak, which in all models has a width of ~ 50 — 100 MeV 
has little chance to be seen directly in this reaction, at least in the considered 
energy region. 
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E°™ = 1 MeV 



40 



60 



E n MeV 
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Figure 5: Effect of decreasing 1^(1405) on the neutron spectra 



5. Conclusions 

Dynamically exact Faddeev-type calculations for the K~ + d — > n + H + n 
reaction were performed in the energy range E^l =0 — 50 MeV in order to 
find the signature of the A(1405) resonance in the observable neutron spectra. 
Four different phenomenological KN — 7rS interactions were used, all well 
reproducing the experimental data in the two-body sector, but having rather 
different pole structure. It was shown, that due to strong kinematical mask- 
ing effect the inclusive neutron spectra do not exhibit a peak, corresponding 
to the A(1405) resonance. We demonstrated, that the deviation spectrum 
method in most of the cases can eliminate the disturbing kinematical factors 
and produce "deviation" spectra clearly showing a resonance peak, which is 
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correlated with the "original" two-body A (1405) resonance. However, the 
shape and position of the peaks in the deviation spectra may significantly 
differ from those of their "original" counterparts. For one of the potentials 
(SIDDHARTA 1) even the deviation spectra show no trace of the resonance, 
probably due to the closeness of the original pole to the KN threshold and 
its large imaginary part. 

Acknowledgment 

The work was supported by the OTKA grant T71989. 

References 

[1] J.Esmaili,Y. Akaishi and T. Yamazaki, Phys. Rev. C 83, 055207, (2011). 

[2] K. Miyagawa and J. Haidenbauer. larXiv: 1202.42721 [nucl-th] (2012). 

[3] N. V. Shevchenko, A. Gal and J. Mares, Phys. Rev. Lett. 98 082301 
(2007), N. V. Shevchenko, A. Gal, J. Mares and J. Revai, Phys. Rev. C 
76, 044004 (2007). 

[4] Y. Ikeda and T. Sato, Phys. Rev. C 76, 035203 (2007). 

[5] N. V. Shevchenko, Phys. Rev. C 85,034001 (2012). 

[6] N. V. Shevchenko. larXiv:1201.3T73 [nucl-th]. (2012) 

[7] A. Bahaoui, C. Fayard, T. Mizutani, and B. Saghai, Phys. Rev. C 68, 
064001 (2003). 

[8] G. Toker, A. Gal, and J. M. Eisenberg, Nucl. Phys. A 362, 405 (1981). 

[9] M.Torres, R.H. Dalitz and A. Deloff, Phys. Lett. B 174,213,(1986) 

[10] T. H. Tan, Phys. Rev. D 7,600,(1973) 

[11] D. Jido, E. Oset and T. Sekihara, Eur. Phys. J. A 42, 257 (2009). 

[12] D. Jido, E. Oset and T. Sekihara, Eur. Phys. J. A 47, 42 (2011). 

[13] A. Alaylioglu, D. Eyre, M. Brannigan and J. P. Svenne, J. Comp. Phys. 
62, 383 (1986). 

16 



